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We consider the possibility that a large fraction of the ultrahigh energy cosmic rays (UHECRs) 
are decay products of Z bosons which were produced in the scattering of ultrahigh energy cosmic 
neutrinos (UHECj/s) on cosmological relic neutrinos (Rz/s). We compare the observed UHECR 
spectrum with the one predicted in the above Z-burst scenario and determine the mass of the 
heaviest Yiu as well as the necessary UHECi/ flux via a maximum likelihood analysis. 



I. INTRODUCTION 

Big-bang cosmology predicts the existence of the cos- 
mic microwave background radiation (CMBR) and a sim- 
ilar background of Ri^s with an average number density 
of {rivi) ~ 56 cm~^ per light neutrino species i {m,y. < 1 
MeV). However, the Rj/s have not been detected until 
now. 

A possibility for their detection was discussed some 
time ago: the UHECz^ spectrum should have absorption 
dips at energies « E'J:^ = M|/(2m^J = 4.2 • lO^^ eV (1 
eV/m^j) due to resonant annihilation with the R^s into 
Z bosons of mass Mz 0H|. Recently it was realized that 
the same annihilation mechanism might already be visi- 
ble in the UHECR spectrum at energies above the 
predicted Greisen-Zatsepin-Kuzmin (GZK) cutoff |^-||] 
around 4-10^^ eV. It was argued that the UHECRs above 
the GZK cutoff are mainly protons from Z decay. 

This hypothesis was discussed in several papers [p|-p^ . 
We report here on our recent quantitative investigation 
of the Z-burst scenario, where we have determined the 
mass of the heaviest Ki/ as well as the necessary UHECi' 
flux via a maximum likelihood analysis jl^] . 

Our comparison of the Z-burst scenario with the ob- 
served UHECR spectrum was done in four steps. First, 
we determined the probability of Z production as a func- 
tion of the distance from Earth. Secondly, we exploited 
collider experiments to derive the energy distribution of 
the produced protons in the lab system. Thirdly, we 
considered the propagation of the protons, i. e. we deter- 
mined their energy losses due to pion and e+e^ produc- 
tion through scattering on the CMBR and due to their 
redshift. The last step was the comparison of the pre- 
dicted and observed spectra and the extraction of the 
mass of the Kv and the necessary UHECj^ flux. 



II. Z-BURST SPECTRUM 



Our prediction of the contribution of protons from Z- 
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bursts to the UHECR spectrum, for degenerate v masses 
(rui, « "rriui), can be summarized as 

oo -Ro oo 

j{E,m,)^I-F^'-Y,JdEpJdrJde (1) 
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F,^ {E,^ , r)n,, {r)a{e)Q{Ep) {-dP{r, Ep, E)/dE) , 

where the total time and angle integrated detector area 
/ and the normalization factor Fz , which is proportional 
to the sum of the ly fluxes at centre-of-mass (CM) en- 
ergy Mz, are determined later by the comparison with 
the UHECR data. E is the energy of the protons ar- 
riving at Earth. Further important ingredients in our 
prediction (0) arc: the UHECi^ fluxes Fi,. {E^. , r) at the 
resonant energy E^. « E^J'^ and at distance r to Earth, 
the number density n^. (r) of the Ri/s, the Z production 
cross section cr(e) at CM energy e — 2 triu E,^. , the en- 
ergy distribution Q{Ep) of the produced protons with 
energy Ep, and the probability P(r, Ep, E) that a proton 
created at a distance r with energy Ep arrives at Earth 
above the threshold energy E. 

The last three building blocks, a, Q, and P, are very 
well determined. At LEP and SLC millions of Z bosons 
were produced and their decays analyzed with extreme 
high accuracy. We used existing published ||l^-[l9[] and 
some improved unpublished pO| data to determine the 
proton momentum distribution Q{Ep). Due to the large 
statistics, the uncertainties of our analysis related to Z 
decay turned out to be negligible. Similarly, the CMBR 
is known to a high accuracy. It plays the key role in 
the determination of the probability P{r, Ep, E) |^ , p2| , 
which takes into account the fact that protons of extra- 
galactic (EG) origin and energies above « 4 ■ 10^® eV lose 
a large fraction of their energies |^,^ due to pion and 
e+e~ production through scattering on the CMBR and 
due to their redshift. P{r, Ep, E), in the form as it has 
been calculated for a wide range of parameters by two of 
the present authors , was an indispensible tool in our 
quantitative analysis. 

Less accurately known in Eq. (|l|) are the first two ingre- 
dients, the flux of UHEC^/s, F,y. {E^. , r), and the neutrino 
number density rij,. (r) . 



The former was assumed to have the form Fi,. (i?^. , r) = 
Fui{E^.,0) (1 + z)", where z is the redshift and where a 
characterizes the source evolution (see also The 
flux at Earth, F^. {E^^ , 0) , has been determined by the 
fit to the UHECR data. In our analysis we went up to 
distances Rq (cf. (0)) corresponding to redshift z — 2 
(cf. [^), and uncertainties of the expansion rate |2j] 
were included. 

The neutrino number density n^. has been treated in 
the following way. For distances below 100 Mpc we varied 
the shape of the n^,. (r) distribution between the homoge- 
neous case and that of mtot(?'), the total mass distribu- 
tion obtained from peculiar velocity measurements [ p5[ . 
In this way we took into account that the density dis- 
tribution of Rfs as hot dark matter (DM) follows the 
total mass distribution; however, with less clustering. It 
should be noted that for distances below 100 Mpc the 
peculiar velocity measurements [ p5[ suggest relative over- 
densities of at most a factor 2 3, depending on the grid 
spacing. We did not follow the unnatural assumption of 
having a relative overdensity of 10^ 10^ in our neigh- 
bourhood, as it was assumed in earlier investigations of 
the Z-burst hypothesis [^|J^,|9|, [lO| , p^ . Our quantitative re- 
sults turned out to be rather insensitive to the variations 
of the overdensities within the considered range, whose 
effect is included in our final error bars. For scales larger 
than 100 Mpc the density was taken according to the 
big-bang cosmology prediction, n^. = 56 ■ (1 -f z)^ cm""^. 



III. DETERMINATION OF M„ 
FLUX 



AND UHECi/ 



We compared the predicted spectrum (|l|) of pro- 
tons from Z-bursts with the observed UHECR spectrum 
(cf. Fig. R). Our analysis included UHECR data of 
AG AS A Fly's Eye |2|-|3^, Havcrah Park |l|,||, 

and HIRES . Due to normalization difficulties we did 
not use the Yakutsk results. 

The predicted number of UHECR events in a bin was 
taken as 



N{^) 



dE [A-E-^ + Fz-jiE,m,)] 



(2) 



where Ei is the lower bound of the i^^ energy bin. The 
first term is the usual power-law behavior, which de- 
scribes the data well for smaller energies [^,^ . For this 
term we studied two possibilities. In the first case we 
assumed that the power part is produced in our galaxy. 
Thus no GZK effect was included for it ("halo"). In the 
second - in some sense more realistic - case we assumed 
that the protons come from uniformly distributed, EG 
sources and suffer from the GZK cutoff ("EG"). In this 
case the simple power-law-like term is modified, by tak- 
ing into account the probability P{r, Ep, E), and falls 
off around 4 • 10^^ eV (see Fig. |^). The second term of 
the fiux in Eq. (||) corresponds to the spectrum of the 
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FIG. 1. The available UHECR data with their error bars 
and the best fits from Z-bursts Note that there are no 

events above 3 x lO'^" eV (shown by an arrow). Bins with no 
events show the 1 a upper bounds on the flux. Therefore the 
experimental value of the integrated flux is in the "hatched" 
region with 68% confldence level ("hatching" is a set of in- 
dividual error bars; though most of them are too large to 
be depicted in full). The dotted line shows the best fit for 
the "halo"-case. The bump around 4 • 10^^ eV is due to the 
Z-burst protons, whereas the almost horizontal contribution is 
the first, power-law-like term of Eq. (^). The solid line shows 
the "extragalactic"-case. The flrst bump at 4- 10^^ eV repre- 
sents protons produced at high energies and accumulated just 
above the GZK cutoff due to their energy losses. The bump 
at 3 • 10^^ eV is a remnant of the Z-burst energy. The dashed 
line shows the contribution of the power-law-like spectrum 
with the GZK effect included. The predicted fall-off for this 
term around 4 • 10^'' eV can be observed. 



Z-bursts, Eq. (Q). A and Fz are normalization factors. 
Note that the following implicit assumptions have been 
made through the form of Eq. (||): i) We have assumed 
(and later checked) that the UHE photons from Z-bursts 
can be neglected, ii) We have assumed that there are 
no significant additional primary UHE proton fluxes be- 
yond the extrapolation of the above power-law. This con- 
straint will be relaxed in a future publication ||3^ ]. 

The expectation value for the number of events in a 
bin is given by Eq. (^. To determine the most probable 
value for we used the maximum likelihood method 
and minimized Q the x^{0i^iEz,m^), 

26.0 

X^= Yl 2[iV(z)-iVo(*) + iVo(z)ln(7Vo(*)/iV(*))], (3) 

i=18.5 

where No{i) is the total number of observed events in the 
i"^ bin. As usual, we divided each logarithmic unit into 
ten bins. Since the Z-burst scenario results in a quite 
small flux for lower energies, the "ankle" is used as a 
lower end for the UHECR spectrum: log(£'niin/eV) ~ 
18.5. Our results are insensitive to the definition of the 
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upper end (the flux is extremely small there) for which 
we choose log(£'max/eV) = 26. The uncertainties of the 
measured energies are about 30% which is one bin. Using 
a Monte-Carlo method we included this uncertainty in 
the final error estimates. 

In our fitting procedure we had four parameters: 
/3, A, Fz and m^. The minimum of the x^iP^ ^z, m^) 
function is Xmin "^i^min which is the most proba- 
ble value for the mass, whereas '^i') = 
Xo(m,/) = Xmin + 1 gives the la (68%) confidence in- 
terval for m^. Here /3', A' , F'^ are defined in such a way 
that the x^(/3j "^i-) function is minimized in (3, A 
and Fz at fixed nii,. 

Qualitatively, our analysis can be understood in the 
following way. In the Z-burst scenario a small Rz^ mass 
needs large E^^^ in order to produce a Z. Large i?™'' re- 
sults in a large Lorentz boost, thus large Ep. In this 
way the shape of the detected energy (E) spectrum de- 
termines the mass of the R;^. The sum of the necessary 
UHECi^ fluxes was then determined from the obtained 
normalization Fz- 

Our best fits to the observed data can be seen in Fig. |^, 
for evolution parameter a = 1. We found a neutrino mass 

oi2Mtl'llfi'lt\ eV for the "halo"- and 0.26^1^1^0^1] 
for the "EG"-case, respectively. The first numbers are 
the la, the numbers in the brackets are the 2a errors. 
This gives an absolute lower bound on the mass of the 
heaviest v of 0.06 eV at the 95% CL. The fits are rather 
good; for 21 non-vanishing bins and 4 fitted parameters 
they can be as low as = 15.1. We determined for 
a wide range of cosmological source evolution (a = 3) 
and Hubble parameter {Hq = 0.64 0.78 km/sec/Mpc) 
and observed only a moderate dependence on them. The 
results remain within the above error bars. We performed 
a Monte-Carlo analysis studying higher statistics. In the 
near future. Auger [ p7| , p8| will provide a ten times higher 
statistics, which reduces the error bars in the neutrino 
mass to « one third of their present values. 

The necessary UHECi^ flux at El°^ has been obtained 
from our fit values of the normalization Fz- We have 
summarized them in Fig. ^, together with some exist- 
ing upper limits and projected sensitivities of present, 
near future and future observational projects. It is ap- 
parent that the flux determination depends much more 
on the evolution uncertainties than the mass determi- 
nation. The necessary i' flux appears to be well below 
present upper limits and is within the expected sensitiv- 
ity of AMANDA, Auger, and OWL. 
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FIG. 2. Differential neutrino fluxes required by the 
Z-burst hypothesis for the "halo" and the "extragalactic" 
case jl^. The evolution parameter a takes values 0,1,2,3 from 
top to bottom for both cases, respectively. The horizontal er- 
rors indicate the 1 a uncertainty of the mass determination 
and the vertical errors include also the uncertainty of the 
Hubble expansion rate. Also shown are upper limits from 
Fly's Eye and the Goldstone lunar ultrahigh energy neu- 
trino experiment GLUE |4Cl||, as well as projected sensitivities 
of AMANDA @, Auger^Q and OWL 



tide. Using the mass difference of the atmospheric v os- 
cillation for the heaviest mass (24j, one obtains values 
between 0.03 and 0.09 eV. It is an intriguing feature of 
our result that the smaller one of the predicted masses is 
compatible on the « 1.3 cr level with this scenario. 

Another popular possibility is to have 4 neutrino types. 
Two of them - electron and sterile neutrinos - are sepa- 
rated by the solar j/ oscillation solution, the other two - 
muon and tau - by the atmospheric v oscillation solution, 
whereas the mass difference between the two groups is of 
the order of 1 eV. We studied this possibility, too. On our 
mass scales and resolution the electron and sterile neu- 
trinos are practically degenerate with mass mi and the 
muon and tau neutrinos are also degenerate with mass 
7712 • The best fit and the 1 a region in the ttii — 7772 plane 
is shown in Fig. |^ for the "EG" -case. The dependence 
of this result on the cosmological evolution and on the 
UHECi^ spectrum will be discussed elsewhere Since 
this two-mass scenario has much less constraints the al- 
lowed region for the masses is larger than in the one-mass 
scenario. 




Auger ^- 



OWL 



uL 



uL 



uL 



uL 



uL 



10'^ 10'* 10^° 10^' 10^ 



IV. COMPARISON WITH AM^ FROM u 
OSCILLATIONS 

One of the most attractive patterns for v masses is 
similar to the one of the charged leptons or quarks: the 
masses are hierarchical, thus the mass difference between 
the families is approximately the mass of the heavier par- 



V. DIFFERENCES WITH RESPECT TO 
YOSHIDA ET AL. (1998) 

Numerical simulations of Z-burst cascades for 777^ 1 
eV, taking into account all known EG propagation ef- 
fects, were performed by Yoshida et al. ||l3|. Based on 
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FIG. 3. The best fit and the 1 a (68% confidence level) 
region in a scenario with two non-degenerate y masses ||l5| . 



necessary UHECz^ flux was found to be consistent with 
present upper limits and detectable in the near future. 
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case studies, relative overdensities of Ri^s ranging from 
20 10^ on a scale of 5 Mpc were argued to be necessary 
in order to get a successful description of the UHECR 
events and rate above the GZK cutoff without violat- 
ing lower energy photon flux limits and without invoking 
inconceivable UHEC;/ fluxes. For such large overdensi- 
ties, most of the UHECRs from Z-bursts originate nearby 
and their attenuation to the Earth can be neglected. In 
our case, with realistic overdensities < 2 3 on scales 
< 100 Mpc, most of the UHECRs from Z-bursts origi- 
nate from cosmological distances. Therefore, despite of 
the fact that by construction the overall rate of UHECRs 
from Z-bursts observed at Earth is the same in both in- 
vestigations, the predicted spectra are quite different. No 
large overdensity is needed to reproduce the data. Note 
that the EG scenario is dominated not by the nearby Z- 
burst but by the pile-up of Z-burst protons due to the 
GZK effect (cf. Fig. |l|). 



VI. CONCLUSIONS 

We reported on a comparison of the predicted spec- 
trum of the Z-burst hypothesis with the observed UHE- 
CR spectrum [|l^. The mass of the heaviest Tiv turned 
out to be = 2.34j:JJ9 eV for halo and 0.26l^ i^ eV for 
EG scenarios. The second mass, with a lower bound of 
0.06 eV on the 95% CL, is compatible with a hierarchical 
V mass scenario with the largest mass suggested by the 
atmospheric v oscillation. The above u masses are in the 
range which can be explored by future laboratory exper- 
iments like the (3 decay endpoint spectrum and the v less 
/3/3 decay [^,0 . They compare also favourably with the 
tau(?) neutrino mass range 0.04 4.4 eV found recently 
from a detailed analysis of the latest CMBR measure- 
ments We analysed a possible two-mass scenario 
and gave the corresponding confldence level region. The 
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